Adaptive compensation systems for mitigating distortion due to nonlinear power amplifiers

ABSTRACT

A method for pre-processing a signal prior to receipt of the signal by a non-linear device (NLD) to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD. The method includes generating a modified amplitude r k  based on an amplitude x k  derived from a digitized envelope R i,k  of the signal, generating a scale factor c k  based on the modified amplitude r k  and the amplitude x k , generating an AM-AM compensated signal based on the scale factor c k , generating a first AM-PM phase compensation value based on the modified amplitude r k , and generating an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value.

STATEMENT REGARDING FEDERALLY SPONSORED RESEARCH OR DEVELOPMENT

This invention was made with government support under Contract No. FA8802-09-C-0001 awarded by the Department of the Air Force. The government has certain rights in the invention.

BACKGROUND

The power amplifier is one of the most important subsystems of modern communication systems. In the case of satellite downlinks, the power efficiency of high-power amplifiers (HPAs) is important because power requirements directly impact size, weight, and cost of the satellite payload. Because the power efficiency is relatively high when the amplifier is operated near the saturation region, as may be the case for satellite links, the HPA frequently must be operated in such a region. However, the nonlinearity of the amplifier in the near saturation region may introduce considerable distortion in the signal to be amplified. As the output power back off is reduced, the signal to distortion power ratio at the amplifier output is correspondingly reduced. This may place a significant restriction on the amount of back off that needs to be introduced, resulting in a loss of the available output power and, equally important, a reduced power conversion efficiency. This results in an increased demand on D.C. power which, in cases where power is supplied by solar panels, may negatively impact the size, weight, and cost of the satellite payload. Schemes by which the amplifier output back off may be minimized while at the same time mitigating the distortion effects of the amplifier may therefore be desirable. Such mitigation becomes even more important when the amplifier input is comprised of relatively high power signals along with relatively low power signals, as in this case distortion due to the strong signals may result in a very poor signal to distortion power ratio for the weak signals.

Previous techniques to mitigate the amplifier distortion effects generally fall into two categories. In the first category may be techniques that attempt to cancel the nonlinear term arising due to the power amplifier in a feedforward or feedback mode with the compensation operation at radio frequency (RF). The feedforward linearizer consists of two loops. The first loop subtracts samples of the input signal from samples of the amplifier output signal to produce samples of the main amplifier's distortion. The second loop subtracts the amplified sampled distortion from the delayed version of the main amplifier output to obtain the final linearized output. This arrangement is complex to implement, requires a second amplifier that needs to be linear to avoid generating its own distortion terms, and results in power loss due to signal combining at the amplifier output. In the feedback linearizer, the amplifier's input and output are detected and low pass filtered and the resulting baseband signals are compared. The error signal is used to control the gain of the amplifier so as to minimize the distortion. This technique suffers from the bandwidth limitation on the amplifier input signal, as the feedback system may respond to frequencies that are much smaller than the inverse of the delay introduced by the amplifier and associated feedback circuitry; thus, the technique is limited to relatively narrowband signals.

In the predistortion linearization technique, the amplifier input is predistorted such that the overall distortion due to the linearizer and amplifier is minimized. The linearizer gain and phase is obtained iteratively for different input power levels. A digital signal processor (DSP) version of the Cartesian predistortion scheme has been proposed wherein a look up table storing inphase and quadrature components of the linearizer as a function of the input signal envelope is used for predistortion for a set of input signal envelope values. The signal to be amplified is digitized and sampled values of the signal are modified by interpolated value of the stored inphase and quadrature components of the linearizer. The correction is limited by interpolation errors. In a proposed adaptive version of the DSP Cartesian predistortion scheme, the stored values are updated according to the amplifier output signal, which is also limited due to the interpolation errors. Based on the feedback architecture of this scheme, a power amplifier linearizer for a time division duplex system has been proposed wherein a receiving subsystem is shared between the receiver and the power amplifier feedback subsystem, resulting in some reduction of complexity for the time division duplex system.

SUMMARY

Methods for pre-processing a signal prior to receipt of the signal by a non-linear device (NLD) to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD are disclosed. In one embodiment, the method includes: generating a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generating a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generating an AM-AM compensated signal based on the scale factor c_(k); generating a first AM-PM phase compensation value based on the modified amplitude r_(k); and generating an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value. In another embodiment, the method includes: compensating, by at least one processor, an AM-AM distortion of the signal based on a scale factor rk/xk·x_(k) is an amplitude derived from a digitized envelope R_(i,k) of the signal, and r_(k) is a modified amplitude determined by a function B(x_(k)). B is a polynomial function of degree N_(a) and an inverse function of an amplitude transfer function h of the NLD. The method further includes compensating, by at least one processor, an AM-PM distortion of the signal based on a first AM-PM phase compensation value determined by a function φ(r_(k)). φ is a polynomial of degree n_(p) approximating the amplitude to phase transfer function of the NLD about a nominal operating amplitude A_(0,n) of r_(k).

Also disclosed is a system that, in one embodiment, includes a non-linear electronic signal processing device for processing an input signal, and a pre-processor in communication with the non-linear electronic signal processing device. The non-linear electronic signal processing device has a non-linear transfer function. The pre-processor receives a signal, processes the received signal, and outputs the processed signal to an input of the non-linear electronic signal processing device. The pre-processor includes at least one processor programmed to: generate a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generate a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generate an AM-AM compensated signal based on the scale factor c_(k), to generate a first AM-PM phase compensation value based on the modified amplitude r_(k); and generate an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value.

DESCRIPTION OF THE FIGURES

Various embodiments of the present invention are described here by way of example in conjunction with the following figures, wherein:

FIGS. 1, 7 and 8 illustrate adaptive compensation systems according to various embodiments;

FIG. 2 illustrates a block diagram of a digital baseband to IF converter according to one embodiment;

FIG. 3 illustrates a block diagram of a nominal operating point estimator according to one embodiment;

FIG. 4 illustrates a block diagram of an AM-PM compensation synthesizer according to one embodiment;

FIG. 5 illustrates a block diagram of an AM-AM compensation synthesizer according to one embodiment;

FIG. 6 illustrates a block diagram of a digital IF to complex baseband converter according to one embodiment;

FIG. 9 illustrates a block diagram of an adaptive phase estimator according to one embodiment;

FIG. 10 is a plot of amplitude to phase transfer characteristics of TWTA and its polynomial approximation;

FIGS. 11 and 12 are plots of the first and second order derivatives, respectively, of a polynomial approximation g(r) of the amplitude to phase characteristic;

FIG. 13 is a plot of an approximation of the φ(r) function with a degree 3 polynomial;

FIG. 14 is a plot of an approximation of the φ(r) function with a degree 2 polynomial;

FIG. 15 illustrates plots of an amplitude transfer characteristics for a Saleh model and an amplitude transfer characteristic compensated by a degree 7 polynomial;

FIG. 16 illustrates plots of an amplitude transfer characteristics for a Saleh model and an amplitude transfer characteristic compensated by a degree 9 polynomial;

FIG. 17 is a plot of amplifier output back off as a function of RF input signal power with the sum of two OQPSK signals at the amplifier input;

FIG. 18 illustrates plots of probability of bit error P_(e) for the WB signal versus (E_(b)/N₀)_(WB) with and without compensation for an output back off of 3 dB;

FIG. 19 illustrates plots of probability of bit error P_(e) for the NB signal versus (E_(b)/N₀)_(NB) with and without compensation for an output back off of 3 dB;

FIG. 20 illustrates plots of probability of bit error P_(e) for the WB signal versus (E_(b)/N₀)_(WB) with and without compensation for an output back off of 2 dB;

FIG. 21 illustrates plots of probability of bit error P_(e) for the NB signal versus (E_(b)/N₀)_(NB) with and without compensation for an output back off of 2 dB;

FIG. 22 illustrates plots of probability of bit error P_(e) for the WB signal versus (E_(b)/N₀)_(WB) with and without compensation for an output back off of 1 dB;

FIG. 23 illustrates plots of probability of bit error P_(e) for the NB signal versus (E_(b)/N₀)_(NB) with and without compensation for an output back off of 1 dB;

FIG. 24 illustrates plots of probability of bit error P_(e) for the WB signal versus (E_(b)/N₀)_(WB) with and without compensation for an output back off of 0.5 dB;

FIG. 25 illustrates plots of probability of bit error P_(e) for the NB signal versus (E_(b)/N₀)_(NB) with and without compensation for an output back off of 0.5 dB; and

FIG. 26 illustrates a computing device according to one embodiment.

DESCRIPTION

Embodiments of the present application are directed to a number of polar predistortion-based linearization architectures wherein the gain and phase may be individually compensated. In many cases, phase distortion is more important than gain distortion and may be substantially or completely eliminated by architectures presented herein. In certain architectures presented herein, the linearizer gain function may be represented by a polynomial of odd degree having coefficients obtained by an explicit minimization of the mean squared error between the actual amplifier output gain response and the ideal response, which may be selected to be linear. The method may also be applied to other selected responses. The minimization may be performed using a gradient algorithm with a power 4 optimization index. Such an approach provides very good response up to the input signal amplitude where the amplifier output is saturated. However, beyond this value of the saturation input, the polynomial-based linearizer may result in an oscillatory response at the amplifier output. To avoid this problem, a modified linearizer may be used wherein the output of the linearizer is clipped to a value equal to the saturation input for any value exceeding the saturation input. This may result in ideal-limiter characteristics for the linearized amplifier. Simulation examples show that a 7^(th) or 9^(th) degree polynomial provides excellent linearization results. Because the amplifier gain response is selected to be an odd function, only the odd degree terms of the linearizer are non zero. Thus, only 4 or 5 coefficients may need to be stored for the linearization compared to a large array of gain values as suggested by previous methods.

According to various embodiments, the amplitude to phase transfer characteristic of the amplifier may be modeled by a polynomial of the input signal envelope. Typically, a 9^(th) or lower degree polynomial may be needed to precisely model the characteristics over the entire range of the input signal envelope. Because the range of input signal may be limited for any given power level, the approximation may be performed by a 2^(nd) or 3^(rd) degree polynomial about a nominal operating point, resulting in a significant computational advantage which may be important in applications involving very wide band signals which may need to be sampled at relatively high rates. The coefficients of the relatively low degree polynomial may be updated at time intervals over which the nominal input power level of the amplifier is expected to change significantly due to various factors such as channel variations. The architecture may include a method to evaluate and update the nominal power level as well. At every sample time, the expected phase distortion due to the amplifier may be evaluated using a 2^(nd) or 3^(rd) degree polynomial and pre-compensated for in the AM-PM compensator. Architectures that operate at RF and baseband are presented. Open loop adaptive architectures may adapt with respect to the nominal input power level changes.

In addition, close loop architectures are presented that may also adapt to the (relatively slow) changes in the amplifier characteristics as well. In these architectures, the error signal between the output of the amplifier with the open loop compensation and the ideal response may be modeled by a parametric model such as a polynomial having coefficients adapted to the variations that may occur in the amplifier characteristics and any other component variations. Because these variations may occur over a much larger time scale compared to the time scale of the signal variations, the error model parameters may remain essentially constant. Therefore, generation of the adaptive close loop compensation signal to offset the variations in the characteristics of the amplifier and any other components need not depend upon the instantaneous phase or amplitude error. The adaptive close loop compensation signal may instead be derived from the amplifier input signal amplitude and phase without using instantaneous feedback from the amplifier. Accordingly, the close loop adaptive architecture of the present application does not have a bandwidth limitation that arises in other feedback architectures of the prior art that require instantaneous output feedback.

Also, in contrast to the Cartesian predistortion architecture of the prior art, in the close loop architecture of the present application there may be two distinct compensation signals based on the known amplifier characteristics and the variations in the latter, respectively. Because the second component of the compensation is relatively insignificant compared to the first, the architecture disclosed in the present application may be relatively more robust in terms of stability compared to architectures having a single compensation signal.

Simulations have shown that embodiments of the adaptive compensation system of the present application provide significant performance improvement, even when applied to the AM-AM and AM-PM characteristics of an otherwise state of the art linearized amplifier. It is believed that these embodiments provide a distortion mitigation of a much higher level than is currently available.

In one general aspect, embodiments of the present application are directed to devices and methods for pre-processing a signal prior to receipt of the signal by a non-linear device (NLD) to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD. AM-AM distortion is caused by a change in a gain of the NLD arising from a change in an input power level of the NLD. AM-PM distortion is caused by a phase shift introduced by the NLD arising from a change in an input power level of the NLD. An NLD that receives an input bandpass signal may comprise only the power amplifier or multiple linear or nonlinear devices such as converters, amplifiers, etc. in a cascaded arrangement. The input bandpass signal may be, for example, a single digitally modulated signal such as an offset quadrature phase-shift keying (OQPSK) signal or the sum of a number of digitally modulated signals. In either case, the combined signal at the input of the NLD may be represented in the following equivalent forms:

The AM-PM form: s _(i)(t)=R _(i)(t)cos(2πf ₀ t+θ _(i)(t))   (1)

The inphase and quadrature form: s _(i)(t)=u _(i)(t)cos(2πf ₀ t)−v _(i)(t) sin(2πf ₀ t)   (2)

where R_(i)(t) and θ_(i)(t) are respectively the real envelope and the phase of the bandpass signal s_(i)(t), and u_(i)(t) and v_(i)(t) respectively denote the real inphase and quadrature baseband signals. The frequency f₀ may be the RF frequency f_(c) when the nonlinear device is just the power amplifier or it may be equal to the intermediate frequency (IF) frequency f_(IF) when the NLD is comprised of both the up converter and RF power amplifier. The signal s_(i)(t) may be available at the complex baseband and is equal to (u_(i)(t)+j v_(i)(t)) with j≡√{square root over (−1)}, in which case the real base band signals u_(i)(t) and v_(i)(t) are available directly as the real and imaginary components of s_(i)(t). As will be appreciated from the discussion that follows, embodiments of the present application may be implemented differently depending on, for example, whether the signal s_(i)(t) is at the RF/IF frequency or at the complex baseband. The output of the NLD, the last stage of which may be the RF power amplifier, may be represented by:

s ₀(t)=R ₀(t)cos(2πf _(c) t+θ _(i)(t)+φ(t))   (3)

where the output real envelope R₀(t) and the phase φ(t) may be obtained in terms of the amplitude transfer characteristics function h( ) and the amplitude to phase transfer characteristics function φ( ). These two functions are directly related to the AM-AM and AM-PM characteristics of the device, which are specified in terms of power levels of the sine wave signal used in the measurements of these characteristics. In terms of the amplitude transfer characteristics function h( ) and the amplitude to phase transfer characteristics function φ( ), the output real envelope R₀(t) and the phase φ(t) may be represented by:

R _(o)(t)=A _(ON) h(x(t));   (4a)

φ(t)=φ(x(t)); x(t)=R _(i)(t)/A _(iN)   (4b)

In equations (4a) and (4b) A_(iN) and A_(oN) may be appropriate input and output normalization amplitudes. For example, A_(iN) and A_(oN) may be selected such that A_(iN)=√{square root over (2P_(i,s))}; A_(oN)=√{square root over (2P_(o,s))} where P_(o,s) denotes the output power at saturation with the sine wave input signal and P_(i,s) denotes the corresponding input power level in watts. However, it will be appreciated that any other convenient normalization may also be used with the corresponding modification in the function h( ) and φ( ). If there are linear devices in the cascade operating at relatively low signal power levels (e.g., bandpass filters), such devices may introduce an amount of frequency dependent phase in the signal which will be added to the phase φ(t) in equation (3). The effect of such devices on the AM-AM and AM-PM model of the composite device is otherwise assumed to be insignificant.

FIG. 1 shows one embodiment of an adaptive compensator 100 that may be used in cases where the inphase and quadrature baseband signals u_(i)(t) and v_(i)(t) are directly available, such as when the compensation is applied at the transmitter or at the baseband processor-type satellite repeaters. The adaptive compensator 100 may be implemented, for example, using a computing device, such as, for example, computing device 2600 discussed below in connection with FIG. 26. As shown in FIG. 1, the signals u_(i)(t) and v_(i)(t) are converted into digital form by analog-to-digital converters (ADCs) 102, 104, respectively. In certain applications, such as those using digital processing at the satellite repeaters, the digital signals u_(i,k) and v_(i,k) may be directly available, in which case no additional ADCs may be required. The real envelope R_(i,k) is obtained from the digital signals u_(i,k) and v_(i,k) at blocks 106, 108 and 110 and may be represented by:

R _(i,k)=√{square root over ([u _(i,k) ² +v _(i,k) ²])}  (5)

The real envelope R_(i,k) is normalized by the selected normalization amplitude A_(iN) in the amplitude normalization block 112 to obtain the normalized amplitude x_(k), which may be represented by:

x _(k) =R _(i,k) /A _(iN) ; k=0, 1, 2, . . .   (6)

The normalized amplitude x_(k) may be input to the AM-AM compensation synthesizer 114, which modifies the input x_(k) by a function B_(m)( ) such that when the modified amplitude r_(k)=B_(m)(x_(k)) is input to the amplitude transfer characteristics function h( ) the output of the amplitude transfer characteristics function h( ) is close to x_(k), i.e., B_(m)( ) achieves amplitude linearization. The AM-AM compensation synthesizer 114 also outputs the scale factor c_(k) which is equal to r_(k)/x_(k) if x_(k)>ε for some relatively small positive number ε and is equal to the derivative of B_(m)(x) at x=0 otherwise, so as to avoid any numerical problems incurred in dividing one nearly zero number by another number also close to zero.

The scale factor c_(k) may be input to the AM-AM compensator 116 where it multiplies the inphase and quadrature baseband signals u_(i,k) and v_(i,k) to generate the AM-AM compensated signals u_(p,k) and v_(p,k), respectively, at the AM-AM compensator 116 output. The modified amplitude r_(k) at the output of the AM-AM compensation synthesizer 114 in FIG. 1 may be input to the AM-PM compensation synthesizer 118 which generates the negative of the estimate of the phase that will be introduced by the NLD, denoted by {circumflex over (φ)}_(c,k). The phase compensation can be input to trigonometric blocks 120, 122 which compute the cosine and sine, respectively, of the phase, denoted by a_(k) and b_(k) respectively, as follows:

a _(k)=cos({circumflex over (φ)}_(c,k)); b _(k)=sin({circumflex over (φ)}_(c,k))   (7)

When the sampling rate is relatively high, the trigonometric functions may be estimated by the following polynomial type approximations wherein the number of terms is dependent upon the maximum absolute value of {circumflex over (φ)}_(c,k). For |{circumflex over (φ)}_(c,k)|≦0.3 rad , the following approximation may be used:

cos(φ)≅1−φ²/2; sin(φ)=φ−φ³/6   (8)

The outputs of the trigonometric blocks 120, 122 (a_(k) and b_(k), respectively) may be input to the AM-PM compensator 124 which also receives the AM-AM compensated inphase and quadrature signals u_(p,k) and v_(p,k). The AM-PM compensator 124 may add the phase {circumflex over (φ)}_(c,k) to the phase of the AM-AM compensated complex baseband signal (u_(p,k)+j v_(p,k)); j≡√{square root over (−1)} according to the following equations:

u _(c,k) =a _(k) u _(p,k) −b _(k) v _(p,k)   (9a)

v _(c,k) =a _(k) v _(p,k) +b _(k) u _(p,k)   (9b)

The AM-AM and AM-PM compensated signals u_(c,k) and v_(c,k) may be converted into analog form using digital-to-analog converters (DACs) 126, 128, respectively, and low pass filtered with the low pass filters (LPFs) 130, 132, respectively. The 3 dB bandwidths of the LPFs 130, 132 may be selected to be somewhat higher than the bandwidth of the input baseband signals u_(i)(t) and v_(i)(t) to take into account the nonlinear compensation introduced by the AM-AM and AM-PM compensators 116, 124, respectively. The filtered baseband signals may be input to the baseband to IF up converter 134, which up converts the complex baseband signal to the IF frequency. After IF filtering and amplification at the IF frequency at block 136, the IF signal may be up converted to the RF frequency at block 138. The resulting RF signal may then be input to the RF power amplifier 140, the output of which is the desired compensated RF signal s_(oc)(t). It will be appreciated that, in certain embodiments, some or all of the blocks following the LPFs 130, 132 that filter the AM-AM and AM-PM compensated signals u_(c,k) and v_(c,k) may not be a part of the adaptive compensator 100 and may instead comprise the NLD for which AM-AM and AM-PM compensation is provided by the adaptive compensator 100.

In another embodiment, the up conversion from the baseband signal to the IF frequency may be performed in the digital domain. Selecting the IF frequency such that f_(IF)=(2κ+0.5)f_(s) for some integer κ≧1 with f_(s) denoting the sampling frequency of the baseband signals u_(i,k) and v_(i,k) in FIG. 1, and delaying the sampling clock (not shown) for the ADC 102 with respect to the sampling clock (not shown) for the ADC 104 in FIG. 1 by T_(s)/2 with the sampling period T_(s)=1/f_(s), the sampled IF signal s_(r)(mT_(sH)) may be represented by:

s _(r)[2mT _(sH)]=(−1)^(m) u _(c,m);

s _(r)[(2m+1)T _(sH)]=−(−1)^(m) v _(c,m) ; m=0, 1, . . .   (10a)

In equation (10a), the sampling interval T_(sH) for the IF signal s_(r) may be equal to 0.5 T_(s) in order to satisfy the Nyquist criteria. The expression in equation (10a) may be derived from equation (10b) below for relating the sampled IF signal to the sampled baseband signals u_(c,k) and v_(c,k) when all of the signals are sampled at a rate equal to 2f_(s).

s _(r)(kT _(sH))=Re{[u _(c,k) +jv _(c,k)]exp[j2πf _(IF) kT _(sH) ]}; j≡√{square root over (−1)}; k=0, 1, . . .   (10b)

In equation (10b), Re denotes the operator that takes the real part of its argument. The equation (10a) is obtained from equation (10b) by substitution for f_(IF) and by evaluating the complex exponential function in equation (10b). FIG. 2 illustrates a block diagram of a digital baseband to IF converter 200 according to one embodiment for implementing such a conversion. Referring to FIG. 2, the digital baseband signals u_(c,k) and v_(c,k) at the output of the adaptive compensator 100 in FIG. 1 are input into a pair of multipliers 202, 204 and are multiplied by the sequence (−1)^(k) output by a sequence generator 206 such that alternate samples of the signals u_(c,k) and v_(c,k) are multiplied by −1. The output of multiplier 202 and the inverted output of multiplier 208 are input to a multiplexer 210 that multiplexes the signals from the multipliers 202, 208 into a single multiplexed signal by alternately selecting samples from its inputs. The multiplexer output signal is input to a DAC 212. The output of the DAC 212 is filtered and amplified by an analog BPF/amplifier 214 to generate the IF signal s_(IF)(t). If the selected IF frequency is equal to f_(IF)=(κ+0.5)f_(s) for some odd integer κ≧1, then the embodiment of FIG. 2 may be modified by deleting the multiplier 208.

Recall that r_(k) denotes the output of the AM-AM compensation synthesizer and is equal to the real envelope of the signal that is input to the IF filter and amplifier. The amplifier amplitude to phase transfer characteristic function may be modeled by a polynomial of degree N_(p) over the entire measurement range. Outside the measurement range the phase contributed by the power amplifier, or by the NLD in general, is assumed to have a constant value. The NLD in FIG. 1 may comprise the IF filter and amplifier 136, upconverter 138 and the RF power amplifier 140, wherein generally only the RF power amplifier 140 introduces the nonlinear distortion. In such cases, the amplitude to phase transfer characteristic of the NLD may thus be the same as that of the RF power amplifier 140. The phase introduced by the RF power amplifier 140 or the NLD at time instance k and denoted by φ_(k) may be represented by:

$\begin{matrix} {\varphi_{k} = \left\{ \begin{matrix} {{g\left( r_{k} \right)};{A_{1} \leq r_{k} \leq A_{h}}} \\ {\varphi_{1};{r_{k} < A_{1}}} \\ {\varphi_{h};{r_{k} > A_{h}}} \end{matrix} \right.} & (11) \end{matrix}$

where A_(l) and A_(h) respectively denote the minimum and maximum value of the input amplitude range over which measurements are available, φ_(l) and φ_(h) denote the phase contribution at these amplitudes, and g(r_(k)) is a polynomial of degree N_(p) that may be represented by:

g(r _(k))=g ₀ +g ₁ r _(k) +g ₂ r _(k) ² + . . . +g _(N) _(p) r_(k) ^(N) ^(p)   (12)

From equation (12) the derivatives of various orders up to some order N_(p) may be evaluated as:

$\begin{matrix} \begin{matrix} {\mspace{79mu} {{{G_{1}(r)} \equiv \frac{{g(r)}}{r}} = {g_{1} + {2\; g_{2}r} + {3\; g_{3}r^{2}} + \ldots + {N_{p}g_{N_{p}}r^{N_{p} - 1}}}}} \\ {= {G_{1,0} + {G_{1,1}r} + {G_{1,2}r^{2}} + \ldots + {G_{1,{N_{p} - 1}}r^{N_{p} - 1}}}} \end{matrix} & \left( {13a} \right) \\ \begin{matrix} {\mspace{79mu} {{{G_{2}(r)} \equiv \frac{^{2}{g(r)}}{r^{2}}} = {{2\; g_{2}} + {6\; g_{3}r} + \ldots + {{N_{p}\left( {N_{p} - 1} \right)}g_{N_{p}}r^{N_{p} - 2}}}}} \\ {= {G_{2,0} + {G_{2,1}r} + {G_{2,2}r^{2}} + \ldots + {G_{2,{N_{p} - 2}}r^{N_{p} - 2}}}} \end{matrix} & \left( {13b} \right) \\ {\mspace{76mu} \vdots} & \; \\ \begin{matrix} {{{G_{j}(r)} \equiv \frac{^{j}{g(r)}}{r^{r}}} = {{\frac{j!}{0!}g_{j}} + {\frac{\left( {j + 1} \right)!}{1!}g_{j + 1}r} + \ldots + {\frac{N_{p}!}{\left( {N_{p} - j} \right)!}g_{N_{p}}r^{N_{p} - j}}}} \\ {= {G_{j,0} + {G_{j,1}r} + {G_{j,2}r^{2}} + \ldots + {G_{j,{N_{p} - j}}r^{N_{p} - j}}}} \end{matrix} & \left( {13c} \right) \end{matrix}$

With G₀(r) equal to g(r), the coefficients of various derivative polynomials may thus be represented by:

$\begin{matrix} {G_{j,i} = \left\{ \begin{matrix} {{\frac{\left( {j + i} \right)!}{i!}g_{j + i}};} & {{0 \leq i \leq {N_{p} - j}};} & {1 \leq j \leq N_{p}} \\ {0;} & {otherwise} & \; \end{matrix} \right.} & (14) \end{matrix}$

The degree N_(p) of the polynomial g(r) may be relatively high to be implemented in real time compensation. To minimize the number of arithmetic operations that need to be performed for real time compensation, a lower degree polynomial approximation may be obtained about the nominal operating point determined on the basis of the average of amplitude r_(k) during any interval of size T_(I) seconds. The interval T_(I) may be selected such that the average envelope power at the input to the NLD does not have any significant variation over this period. FIG. 3 shows one embodiment of a nominal operating point estimator 300 with input r_(k) providing the estimate of the nominal operating point at the end of the n^(th) averaging interval denoted by A_(0,n) at the output. With Â_(a,n) denoting the estimate of the average value of r_(k) over the n^(th) interval of duration T_(I) seconds over which the signal power is essentially constant, Â_(a,n) may be evaluated by the following equations (15)-(18):

S _(k)=λ_(s) S _(k−1) +r _(k) ; S ₀=0; k=1, 2, . . .   (15a)

A _(a,k) =S _(k)(1−λ_(s))/(1−λ_(s) ^(k+1)); k=1, 2, . . .   (15b)

or equivalently:

A _(a,k)=λ_(s)[(1−λ_(s) ^(k))/(1−λ_(s) ^(k+1))]A _(a,k−1)+[(1−λ_(s))/(1−λ_(s) ^(k+1))]r _(k) ; k=1, 2, . . .   (16)

with 0<λ_(s)<1. The parameter λ_(s) in equations (15)-(16) determines the averaging interval, which may be approximately equal to 1/(1−λ_(s)) samples. With the sampling period equal to t_(s) (sampling rate f_(s)=1/t_(s)), the averaging interval in seconds denoted by T_(I) may be represented by:

T_(I)≅t_(s)/(1−λ_(s))   (17a)

Thus for any specified averaging interval, the desired value of λ_(s) may be obtained by:

λ_(s)=(1−t _(s) /T _(I))   (17b)

In the steady state with k>>N_(s)=1/(1−λ_(s)), the equation (16) may be replaced by the following equation:

A _(a,k)=λ_(s) A _(a,k−1)+(1−λ_(s))r _(k)   (18)

The estimates A_(a,k) may be sampled by a rate f_(I)=1/T_(I) sampling clock 302 at the end of N_(s) sampling intervals, providing the estimates Â_(a,n) for n=0, 1, . . . . The estimate Â_(a,n) at the end of the n^(th) averaging interval may thus be represented by:

Â _(a,n) =A _(a,n) ₀ ; n ₀=(n+1)Ns−1; n=0, 1, 2, . . .   (19)

In another embodiment, instead of the exponential data weighted averaging performed by equations (16)-(18), an arithmetic average may be performed over consecutive segments of N_(s) samples, with N_(s) denoting the integer part of (T_(I)/t_(s)), to obtain Â_(a,n) for n=0, 1, 2, . . . . As shown in FIG. 3, along with the average value of r_(k) over the interval of N_(s) samples as estimated by equations (15)-(19), the estimates of the minimum and maximum values of r_(k) over the 0^(th) interval of N_(s) samples, denoted by Â_(m,0) and Â_(M,0), respectively, may be obtained from equations (20a-c) set forth below. Similar estimates for the n^(th) interval of N_(s) samples, denoted by Â_(m,n) and Â_(M,n), are obtained for n=1, 2, . . . .

A _(m,k)=min(A _(m,k−1) , r _(k)); A _(m,0) =r ₀; 0<k<N _(s)   (20a)

A _(M,k)=max(A _(M,k−1) , r _(k)); A _(M,0) =r ₀; 0<k<N _(s)   (20b)

Â _(m,0) =A _(m,N) _(s) ⁻¹ ; Â _(M,0) =A _(M,N) _(s) ⁻¹   (20c)

The embodiment of FIG. 3 shows an implementation of equations (18)-(20), with the output of blocks 304, 306, and 308 corresponding to Â_(a,n), Â_(m,n) and Â_(M,n), respectively. In FIG. 3, A_(∞) denotes some large number which is an upper bound on the value of r_(k). At the end of the N_(s) samples intervals, the nominal mean value of the amplitude r_(k) may determined by the nominal operating point evaluator 310 in FIG. 3 as follows:

$\begin{matrix} {A_{0,n} = \left\{ \begin{matrix} {{\hat{A}}_{a,n};{{\hat{A}}_{m,n} \geq {A_{1}\mspace{14mu} {and}\mspace{14mu} {\hat{A}}_{M,n}} \leq A_{h}}} \\ {{{0.5\left( {A_{1} + {\hat{A}}_{M,n}} \right)};{{\hat{A}}_{m,n} < {A_{1}\mspace{14mu} {and}\mspace{11mu} {\hat{A}}_{M,n}} \leq A_{h}}}\;} \\ {{0.5\left( {{\hat{A}}_{m,n} + A_{h}} \right)};{{\hat{A}}_{m,n} > {A_{1}\mspace{14mu} {and}\mspace{14mu} {\hat{A}}_{M,n}} > A_{h}}} \end{matrix} \right.} & (21) \end{matrix}$

For the purpose of AM-PM compensation, the phase introduced by the amplifier 140 may be modeled by a relatively low degree polynomial of degree n_(p) around the nominal input amplitude A₀ computed from equation (21). In certain embodiments n_(p) may be in the range of 1 to 3 compared to N_(p), which may be of the order of 10, where N_(p) is the degree of a polynomial that models the amplifier amplitude to phase transfer characteristics over the entire range of the amplitudes.

The desired polynomial approximation of the phase due to the amplifier may be represented by:

$\begin{matrix} {{\varphi (r)} = {{p_{0} + {p_{1}\left( {r - A_{0}} \right)} + {p_{2}\left( {r - A_{0}} \right)}^{2} + \ldots + {p_{n_{p}}\left( {r - A_{0}} \right)}^{n_{p}}} \equiv {p\left( {r - A_{0}} \right)}}} & (22) \end{matrix}$

where the suffix n on A has been dropped for convenience. The coefficients of the polynomial p(r) may be computed in the polynomial φ(r) coefficients evaluator 142 from the polynomial G_(j)(r) specified by equations (13)-(14), as in equations (23)-(24) below:

$\begin{matrix} {{{p_{j} = {{\frac{1}{j!}\frac{{g^{j}(r)}}{r^{j}}}_{r = A_{0}}}};{j = 0}},1,\ldots \mspace{14mu},n_{p}} & (23) \\ {{{p_{j} = {\frac{1}{j!}{G_{j}\left( A_{0} \right)}}};{j = 0}},1,\ldots \mspace{14mu},n_{p}} & (24) \end{matrix}$

The polynomial functions G_(j)(A₀) may be evaluated from their coefficients given by equation (14). Denoting by C_(s) the time in terms of the number of samples it takes to compute the coefficients p_(j) from equations (21) and (24), then after the initial (N_(s)+C_(s)) samples, these coefficients may be kept fixed for a subsequent period of N_(s) samples and used to evaluate φ(r) from equation (22) during this period. The processes of averaging given by equations (15)-(20), and those of computing the coefficients p_(j)'s from equations (21) and (24), may be performed concurrently during the consecutive intervals of N_(s) samples.

Referring to FIG. 1, the nominal operating point A₀ estimated by the nominal operating point estimator 300 and the polynomial coefficients p₀, p₁, . . . , p_(n) _(p) evaluated by the polynomial φ(r) coefficients evaluator 142 may be input to the AM-PM compensation synthesizer 118 at the end of each averaging interval of N_(s) samples. These coefficients may be used by the AM-PM compensation synthesizer 118 for a subsequent period of N_(s) samples. Over any of the consecutive interval of N_(s) samples, the phase compensation {circumflex over (φ)}_(c,k)(r_(k)) required to compensate for the amplifier introduced phase may be computed by the AM-PM compensation synthesizer 118.

FIG. 4 shows one embodiment of an AM-PM compensation synthesizer 400. The AM-PM compensation synthesizer 400 may be identical to the AM-PM compensation synthesizer 118 of FIG. 1 in certain embodiments and may determine the phase compensation {circumflex over (φ)}_(c,k)(r_(k)) in accordance with equation (25) below. The polynomial approximation of the phase due to the amplifier is computed at block 402 and, along with A_(l), A_(h), φ_(l), φ_(h) and the outputs of threshold blocks 404, 406, is input to multiplexer block 408. The output of the multiplexer block 408 is multiplied by −1 to generate the phase compensation {circumflex over (φ)}_(c,k)(r_(k)). The constants A_(l), A_(h), φ_(l), and φ_(h) may be determined for the NLD according to equation (11).

$\begin{matrix} {{{\hat{\varphi}}_{c,k}\left( r_{k} \right)} = \left\{ \begin{matrix} {{- {p\left( {r_{k} - A_{0}} \right)}};{A_{1} \leq r_{k} \leq A_{h}}} \\ {{- \varphi_{1}};{r_{k} < A_{1}}} \\ {{- \varphi_{h}};{r_{k} > A_{h}}} \end{matrix} \right.} & (25) \end{matrix}$

In the case of certain linearized amplifiers, the AM-AM is well linearized, however, there may be a significant AM-PM which may significantly degrade the performance of the digital signal at the amplifier output. Thus, only the AM-PM compensation may be adequate. In such cases, the input to the AM-PM compensation synthesizer 400 r_(k) may be set equal to x_(k), and c_(k) at the AM-AM compensation synthesizer 114 may be set equal to 1 in FIG. 1, thus effectively bypassing the AM-AM compensation synthesizer 114.

For applications where both AM-AM and AM-PM are significant, it may be desirable to compensate for both of these distortion effects. Compensation for the AM-AM distortion may be performed by the AM-AM compensation synthesizer 114 and the AM-AM compensator 116 in FIG. 1. Compensation for AM-AM may be implemented using an inverse function such that the cascade of the inverse function and the amplitude transfer characteristic function yields a linear function. Let h(x) denote the amplitude transfer characteristics function of the amplifier 140 so that the output amplitude y_(k)=h(x_(k)) with the input signal amplitude equal to x_(k). Let the inverse function of h( ) be modeled by a polynomial of degree N_(a) given by:

B(x _(k))=B ₀ +B ₁ x _(k) +B ₂ x _(k) ² + . . . +B _(N) _(a) x _(k) ^(N) ^(a)   (26)

The output amplitude of the cascade of the inverse function and the amplifier 140 may thus be represented by:

y _(k) =h(B(x _(k)))   (27)

The difference between y_(k) and x_(k) is the error e_(k) which needs to be minimized in the approximation of the inverse function. Selecting the error measure to be F(e_(k)) for some nonnegative and monotone increasing function F(e_(k)), for example F(e_(k))=e_(k) ², the minimization problem may then involve the minimization of the following optimization function J with respect to the coefficients of the polynomial B(x), e.g.,

$\begin{matrix} {{\min } = {\int_{0}^{A_{L}}{F\left\{ {{h\left( {{B(x)} - x} \right\}}\ {x}} \right.}}} & (28) \end{matrix}$

where A_(L) is appropriately selected from the amplitude transfer characteristic function h( ). For example, A_(L) may be set equal to the minimum (normalized) input amplitude for which the amplifier output amplitude is equal to the saturation value A_(s). To minimize the optimization function J, the partial derivatives of J with respect to B_(j) for j=0, 1, . . . , N_(a) are set equal to 0. The derivative of J with respect to B_(j) may be represented by:

$\begin{matrix} {\frac{\partial}{\partial B_{j}} = {\int_{0}^{A_{L}}{F^{\prime}\left\{ {{h\left( {B(x)} \right)} - x} \right\} \ {h^{\prime}\left( {B(x)} \right)}x^{j}{x}}}} & (29) \end{matrix}$

where F′( ) and h′( ) denote the derivatives of the functions F( ) and h( ) respectively. The derivative of h( ), when not available in closed form, may be approximated by the first-order difference given in equation (30):

$\begin{matrix} {{h^{\prime}(y)} \cong \frac{{h\left( {y + \Delta} \right)} - {h\left( {y - \Delta} \right)}}{2\Delta}} & (30) \end{matrix}$

for some relatively small Δ. The integral in equation (29) may also be approximated by a summation as follows:

$\begin{matrix} {{\frac{\partial}{\partial B_{j}} \cong {\sum\limits_{k = 0}^{K - 1}\; {{4\left\lbrack {{h\left( {B\left( x_{k} \right)} \right)} - x_{k}} \right\rbrack}^{3}{h^{\prime}\left( {B\left( x_{k} \right)} \right)}x_{k}^{j}\delta}}};{x_{k} = {k\; \delta}};{\delta = {A_{L}/K}}} & (31) \end{matrix}$

where in the example considered in equation (31) the function F(e) is equal to e⁴ and K is some sufficiently large integer. The optimum solution for the coefficients of the inverse polynomial B(x) may be obtained by the solution of the following set of equations:

$\begin{matrix} {{{\frac{\partial}{\partial B_{j}} = 0};{j = 0}},1,{\ldots \mspace{14mu} \left( {N_{a} - 1} \right)}} & (32) \end{matrix}$

The set of equations (32) may be solved by the gradient algorithm:

$\begin{matrix} {{{B_{j,k} = {{B_{j,{k - 1}} - {\mu \frac{\partial}{\partial B_{j}}}}_{B_{j,{k - 1}}}}};{k = 1}},2,{\ldots \mspace{14mu};{j = 0}},1,\ldots \mspace{14mu},N_{a}} & (33) \end{matrix}$

In equation (33) B_(j,k) denotes the estimate of B_(j) at the iteration number k, and μ is some relatively small adaptation parameter. The iteration may terminate when the function h(B(x)) provides a suitable linear approximation. The initial estimates for B_(j) coefficients in equation (33) may be taken equal to:

B_(0,0)=1; B_(j,0)=0 for j>0   (34)

While most appropriate selections of the function F(e) will yield suitable results in terms of the linearization of the amplitude response, simulation examples suggest that a desirable result may be obtained with the following selection for the function F(e):

F(e)=e ⁴+0.25e ⁶   (35)

As may be anticipated, while the linearized amplitude response h(B(x)) has good linearity over the range of x below the saturation value, it exhibits ripples above the saturation values. Moreover, the power amplifier response by itself may exhibit some ripple effects in the saturation region. The function B(x) obtained by the optimization of J in equation (28) may thus be modified as follows. Denoting by A_(th) the minimum value of the input amplitude x for which h(B(x)) is equal to or nearly equal to the amplifier (normalized) output saturation amplitude A_(s), the compensated response B_(m)(x) may be modified as follows:

$\begin{matrix} {{B_{m}(x)} = \left\{ \begin{matrix} {{B(x)};{x \leq A_{th}}} \\ {A_{p};{x > A_{th}}} \end{matrix} \right.} & (36) \end{matrix}$

where A_(p) is the minimum value of x for which h(A_(p))≅A_(s). The combined response h(B_(m)(x)) may thus be approximately linear over the range 0<x<A_(th) and be equal to A_(s) for x>A_(th), thus in effect obtaining a limiter amplifier. One advantage of this arrangement is an increase in the amplifier output power while at the same time achieving a higher linearization.

To minimize the number of real-time operations, the degree N_(a) of the polynomial B(x) may be made dependent upon the average input signal A_(a), and thus, dependent upon the amplifier output back off. As the polynomial B(x) may be an odd degree polynomial, the number of non-zero coefficients may be equal to (N_(a)+1)/2. In certain embodiments, the polynomials corresponding to N_(a)=3, 5, 7 and 9 may be evaluated and stored in memory, and, depending upon the average value A_(a) or the output back off level, the appropriate coefficients may be used. Alternatively, the highest degree polynomial may be used for all times at the cost of a possible increase in the number of real time computations.

Referring to FIG. 1, the polynomial coefficients B₀, B₁, . . . , B_(Na) may be computed at block 144 and as described above and input to the AM-AM compensation synthesizer 114 which may compute the modified amplitude r_(k)=B_(m)(x_(k)) based on the polynomial B(x). The AM-AM compensation synthesizer 114 may also output the scale factor c_(k), which is equal to r_(k)/x_(k) if x_(k)>ε for some relatively small positive number ε and otherwise equal to the derivative of B_(m)(x) at x=0 to avoid any numerical problems incurred in dividing one nearly zero number by another number also close to zero. Referring to FIG. 1, the AM-AM compensator 116 may multiply the inphase and quadrature baseband signals u_(i,k) and v_(i,k) by the scale factor c_(k) to generate the AM-AM compensated baseband signals u_(p,k) and v_(p,k). The AM-AM compensated baseband signals u_(p,k) and v_(p,k) may be input to the AM-PM compensator 124 for compensating for the AM-PM distortion.

FIG. 5 shows the block diagram of the AM-AM compensation synthesizer 500 according to one embodiment. The AM-AM compensation synthesizer 500 may be identical to the AM-AM compensation synthesizer 114 of FIG. 1 in certain embodiments. Referring to FIG. 5, the polynomial B(x) may computed for the input x=x_(k) for the example case of N_(a)=7, for example. The order of computation in FIG. 5 requires a total of seven multiplications (block 502) to generate B(x_(k)). In FIG. 5, the amplitude x_(k) may be compared with the threshold A_(th) at threshold block 504. The output of the threshold block 504 V_(H) may be equal to 1 if x_(k)>A_(th) and be equal to 0 otherwise. Multiplexer block 506 may output A_(p) if V_(H)=1 and output B(x_(k)) otherwise. The output r_(k) may thus equal B(x_(k)) if x_(k)≦A_(th) and equal A_(p) otherwise. The AM-AM compensation synthesizer 500 of FIG. 5 may also compare x_(k) with a relatively small positive number ε in threshold block 508, with the output of the threshold block 508 V_(L) equal to 1 if x_(k)>ε and equal to 0 otherwise. The output of switch S₁ 510 may be equal to 1 if V_(L) is equal to 0 and be equal to x_(k) otherwise. The output of switch S₁ 510 may divide x_(k) in the divider 512, the output of which is connected to switch S₂ 514. The switch S₂ 514 output c_(k) may thus equal B₁ if V_(L) is equal to 0 and equal B_(m)(x_(k))/x_(k) otherwise. The logic used in the generation of c_(k) avoids any numerical problems incurred when dividing one nearly zero number by another number also close to zero. Referring to FIG. 1, the output r_(k) of the AM-AM compensation synthesizer 114 may be input to the AM-PM compensation synthesizer 118, while the output c_(k) may be input to the AM-AM compensator 116.

In cases in which the signal s_(i)(t) is at the IF frequency, the signal s_(i)(t) may be sampled directly to provide the digital inphase and quadrature baseband signals u_(i,k) and v_(i,k). This may be achieved by sampling the IF signal at a rate F_(s) that is related to f_(IF) such that the samples of s_(i)(t) relate directly to u_(i,k) and v_(i,k). For example, selection of a sample rate such that (f_(IF)/F_(s))=(2m+1)/4 for some appropriate integer m provides the desired result as shown below.

With

s _(i)(t)=u _(i)(t)cos(2πf _(IF) t)−v _(i)(t)sin(2πf _(IF) t)   (37)

and with f_(IF)=F_(s)(2m+1)/4, T_(F)=1/F_(s), the sample values of the cos and sin terms are given by:

$\begin{matrix} {{\cos \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\cos \left\lbrack {\pi \; {{i\left( {{2\; m} + 1} \right)}/2}} \right\rbrack}} & \left( {38a} \right) \\ {{{{{\sin \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\sin \left\lbrack {\pi \; {{i\left( {{2\; m} + 1} \right)}/2}} \right\rbrack}};{i = 0}},1,{2\ldots}}{with}} & \left( {38b} \right) \\ \begin{matrix} {{{\cos \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\cos \left\lbrack {\pi \; {k\left( {{2\; m} + 1} \right)}} \right\rbrack}};{i = {2\; k}}} \\ {= \left\{ \begin{matrix} {{1;{k\mspace{14mu} {even}};{i = 0}},4,8,\ldots} \\ {{{- 1};{k\mspace{14mu} {odd}};{i = 2}},6,\ldots} \end{matrix} \right.} \end{matrix} & \left( {39a} \right) \\ \begin{matrix} {{{\cos \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\cos \left\lbrack {{\pi \left( {{2\; k} + 1} \right)}{\left( {{2\; m} + 1} \right)/2}} \right\rbrack}};{i = \left( {{2\; k} + 1} \right)}} \\ {= 0} \end{matrix} & \left( {39b} \right) \\ \begin{matrix} {{{\sin \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\sin \left\lbrack {\pi \; {k\left( {{2\; m} + 1} \right)}} \right\rbrack}};{i = {2\; k}}} \\ {= 0} \end{matrix} & \left( {39c} \right) \\ \begin{matrix} {{{\sin \left( {2\pi \; f_{IF}{iT}_{F}} \right)} = {\sin \left\lbrack {{\pi \left( {{2k} + 1} \right)}{\left( {{2\; m} + 1} \right)/2}} \right\rbrack}};{i = {{2\; k} + 1}}} \\ {{= \begin{Bmatrix} {{{- 1};{k\mspace{14mu} {even}};{i = 1}},5,\ldots} \\ {{1;{k\mspace{14mu} {odd}};{i = 3}},7,\ldots} \end{Bmatrix}};{m\mspace{14mu} {odd}}} \end{matrix} & \left( {39d} \right) \end{matrix}$

In equation (39d), the sign of the terms is reversed if m is even. Thus, from equations (37)-(39), s_(i)(iT_(F))=±u_(i)(iT_(F)) for i even and s_(i)(kT_(F))=±(−v_(i)(kT_(F))) for i odd with alternate signs as per equation (39) with the sign of the first term + for both the even and the odd terms subsequences for m odd.

FIG. 6 shows one embodiment of a digital IF to complex baseband converter 600 to obtain the digital inphase and quadrature baseband signals u_(i,k) and v_(i,k). As shown in FIG. 6, the signal s_(i)(t) may be sampled at a rate F_(s) such that (f_(IF)/F_(s))=(2m+1)/4 with m some odd integer. For example, if m=3, then F_(s)=(4/7)f_(IF) and f_(s)=F_(s)/2. The sampling rate f_(s) should also be adequate to satisfy the Nyquist sampling theorem for bandpass signals, e.g., f_(s) should be greater than or equal to the IF bandwidth. In FIG. 6, the input IF signal s_(i)(t) is sampled and converted into digital form by an ADC 602 generating s_(i)(iT_(F)) at the output. The resulting sampled signal may be input to a serial to parallel converter 604, which splits the digital signal into two streams with even numbered samples in the upper branch and the odd numbered samples in the lower branch in FIG. 6. Both the digital streams may be multiplied by a sequence of alternate +1 and −1 (blocks 606 and 608, respectively) generated by an oscillator 610. The lower stream is then multiplied by −1 at block 612 for the case of odd m. The upper stream then constitutes the sequence u_(i,k), while the lower stream is the sequence v_(i,k). The two signals u_(i,k) and v_(i,k) may be input to the adaptive compensator 100 of FIG. 1, in which case the ADCs 102, 104 for converting the analog baseband signals u_(i)(t) and v_(i)(t) into the digital signals u_(i,k) and v_(i,k) are no longer required.

FIG. 7 shows an adaptive compensator 700 according to an alternative embodiment in which the AM-AM and AM-PM compensation may be applied directly at the RF or IF depending upon the specific system configuration and the particular IF and RF frequencies. The adaptive compensator 700 may be implemented, for example, using a computing device, such as, for example, computing device 2600 discussed below in connection with FIG. 26. The arrangement of FIG. 7 may be, for example, that of a bent pipe type satellite. The received RF signal, after filtering and amplification at block 702, may be input to a coupler 704. The coupler 704 may input a small fraction of the power to a down converter 706 and to an IF band pass filter (BPF) 708 for the purpose of synthesizing the compensation signal. The output of the down converter 706 after filtering may be applied to an envelope detector 710. The envelope detector output 710 R(t) may be digitized with ADC 712 to provide the digitized envelope R_(i,k) at the output. In certain embodiments, depending upon the actual frequencies involved, the envelope detector 710 may be implemented by an analog circuit. Alternatively, the envelope detector 710 and ADC 712 may be replaced by a digital IF to complex baseband converter (such as that shown in FIG. 6) to provide the digital complex baseband signals u_(i,k) and v_(i,k), which may be processed according to equation (5) to generate R_(i,k). The digital real envelope signal R_(i,k) may be normalized by the normalization constant A_(iN) at block 714 as in the compensator 700 of FIG. 1 at block 112, except that it will now include the gain of the coupler 704, IF down converter 706 and IF BPF 708, which are the components in the path between the output of the RF BPF/amplifier 702 and the input of the envelope detector 710. The normalized envelope signal x_(k) may be input to the AM-AM compensation synthesizer 500 to provide a modified envelope signal r_(k) and the gain c_(k)=r_(k)/x_(k).

Operation of the AM-AM compensation synthesizer 500 of FIG. 7 is shown in FIG. 5 and may be same as that for the system 100 of FIG. 1. The gain signal may be converted into the analog form by a DAC 716 and filtered by a LPF 718 to generate the analog gain signal c(t). The analog gain signal c(t) may be input to the product modulator 720. The product modulator 720 may also be connected to the output of the RF BPF/amplifier 702 via the coupler 704. The signal c(t) amplitude may modulate the RF signal input in the product modulator 720, thus multiplying the real envelope by c(t). The product modulator 720 should have a good linearity in terms of both of its inputs. The other output r_(k) of the AM-AM compensation synthesizer 500 is input to the AM-PM compensation synthesizer 400 that generates the phase compensation signal {circumflex over (φ)}_(c,k). The operation of the AM-PM compensation synthesizer 400 of FIG. 7 is shown in FIG. 4 and may be the same as that of the AM-PM compensation synthesizer 118 in the system 100 of FIG. 1. The phase compensation signal may be converted into the analog form by a DAC 722 and filtered by a LPF 724 to generate the analog phase compensation signal {circumflex over (φ)}_(c)(t). The signal {circumflex over (φ)}_(c)(t) may be input to the phase modulator 726 where it modifies the phase of the RF signal at its input by {circumflex over (φ)}_(c)(t). Thus, with s_(i)(t) denoting the RF input to the product modulator 720 and given by equation (40):

s _(i)(t)=R _(i)(t)cos(2πf _(c) t+θ _(i)(t))   (40)

where f_(c) denotes the carrier frequency of the RF signal, the output of the phase modulator 726 s_(ic)(t) may be represented by:

s _(ic)(t)=c(t)R _(i)(t)cos(2πf _(c) t+θ _(i)(t)+{circumflex over (φ)}_(c)(t))   (41)

The compensated signal s_(ic)(t) is input to the power amplifier chain 728 (which may comprise power amplifier 140) to generate the compensated and amplified signal s_(oc)(t) at the output.

The adaptive compensators 100, 700 of FIGS. 1 and 7 adapt the compensator with respect to the power level variations at the input of the NLD. These variations are measured in terms of the nominal operating point A_(0,n), which is the measure of average power over some selected time segment. However, there may also be variations of the AM-AM and AM-PM characteristics of the NLD that may be caused by variations in the environmental temperature and with aging. In the case of a power amplifier, the variations in AM-AM characteristics may be relatively insignificant while the variations in AM-PM characteristics may be somewhat significant. For a more precise compensation, it may be desirable to take into account such variations. FIG. 8 shows one embodiment of closed loop adaptive compensator 800 that takes such variations into account. The adaptive compensator 800 may be implemented, for example, using a computing device, such as, for example, computing device 2600 discussed below in connection with FIG. 26. The compensator 800 of FIG. 8 may implement adaptation with respect to the variations in the AM-PM characteristics. If required, the adaptation with respect to the variations in the AM-AM characteristics can be similarly included in system FIG. 8.

FIG. 8 shows one embodiment of a closed loop adaptive compensator 800 for the mitigation of the distortion effects of the NLD. One difference between the open loop system of FIG. 7 and the closed loop system of FIG. 8 is the introduction of the adaptive feedback subsystem 802 to the system of FIG. 7. The adaptive feedback subsystem 802 of FIG. 8 measures the difference between the phase of the reference signal s_(IF)(t) and the phase of the signal at the NLD output. The information may be used to offset such a phase difference by a correction signal applied to the amplifier input signal. This correction may be in addition to the phase compensation applied by the open loop system. As shown in FIG. 8, the signal s_(IF)(t), after delay equalization at block 804 to account for the delay incurred in the evaluation of the phase correction signal and the delay introduced by the NLD, may be input to a band pass limiter 806 to eliminate any amplitude variations in the signal. The output of the band pass limiter 806 denoted by s₁(t) may thus be represented by:

s ₁(t)=κ₁ cos(2πf _(IF) t+θ _(i)(t))   (42)

for some constant κ₁. The output of the power amplifier chain 728 may be sensed via a coupler 808 and applied to a down converter 810 to shift the spectrum to IF. After filtering by a BPF 812, the output of the down converter 810 may be shifted in phase by π/2 at block 814. The phase shifted signal may be input to a second band pass limiter 816 having an output s₂(t) represented by:

s ₂(t)=κ₂ cos(2πf _(IF) t+θ _(i)(t)+ψ(t)−{circumflex over (ψ)}(t))   (43)

for some constant κ₂. In equation (43), ψ(t) may denote the phase error incurred in the open loop adaptation due to variations in temperature etc., and {circumflex over (ψ)}(t) may be the feedback due to the adaptive feedback subsystem 802. The output of the second band pass limiter 816 s₂(t) may be input to the product demodulator 818, along with the output of the band pass limiter 806. The output of the product demodulator 818 may be represented by:

ξ(t)=κ₀ sin(ψ(t)−{circumflex over (ψ)}(t))   (44)

for some constant κ₀. The signal ξ(t) may be digitized by an ADC 820 to generate ξ_(k−1). The phase error (ψ_(k−1)−{circumflex over (ψ)}_(k−1)) is expected to be relatively small due to the open loop correction already incorporated, and thus the sampled error signal ξ_(k−1) at time (k−1) may be approximated by:

ξ_(k−1)=κ₀(ψ_(k−1)−{circumflex over (ψ)}_(k−1))   (45)

The signal ξ_(k−1) and r_(k) may be input to an adaptive phase estimator 822 to generate the phase estimate {circumflex over (ψ)}_(k), which may then be subtracted from the phase compensation signal {circumflex over (φ)}_(o,k) generated by the open loop compensator 700 at block 824 to create the resultant phase compensation signal {circumflex over (φ)}_(c,k). The resultant phase compensation signal, after being converted into the analog form by the DAC 722 and filtered by the LPF 724, may be applied to the phase modulator 726 to compensate for the phase distortion introduced by the NLD.

FIG. 9 shows a block diagram of an adaptive phase estimator 900 according to one embodiment. In the adaptive phase estimator 900, the phase distortion ψ_(k) arising due to the variations in the AM-PM characteristics of the NLD or due to any other unmodeled non-linearity and/or source of distortion may be estimated by a polynomial q( ) of the amplitude r_(k), wherein the degree of the polynomial q( ) denoted n_(q) is relatively small. The estimate of ψ_(k) denoted by {circumflex over (ψ)}_(k) may thus be represented by:

{circumflex over (ψ)}_(k) =q ₀ +q ₁ r _(k) +q ₂ r _(k) ² + . . . +q _(n) _(q) r _(k) ^(n) ^(q)   (46)

The polynomial coefficients q₀, q₁, . . . , q_(nq) may be selected so as to minimize some measure of the error ε_(k)=ψ_(k)−{circumflex over (ψ)}_(k). Selecting a quadratic performance measure, the polynomial coefficients may be adaptively estimated using the stochastic gradient algorithm. Thus,

{circumflex over (q)} _(j,k) ={circumflex over (q)} _(j,k−1) +μκr _(k−1) ^(j)ε_(k−1) ; j=0, 1, . . . , n _(q) ; k=0, 1, 2 . . .   (47)

In FIG. 9, the polynomial coefficients are estimated at block 902, and {circumflex over (ψ)}_(k) is determined at block 904. The initial estimate of the polynomial coefficients may be all set to 0. FIG. 9 shows the adaptive phase estimator 900 described by equations (46) and (47) for the example case of n_(q)=2. As may be inferred from FIG. 9, the inputs to the adaptive phase estimator 900 are ξ_(k−1)=κ(ψ_(k−1)−{circumflex over (ψ)}_(k−1)) and the amplitude signal r_(k), with {circumflex over (ψ)}_(k) generated at the output. In certain embodiments, only the closed loop AM-PM adaptation may be implemented without the open loop component of the AM-PM adaptation.

The performance of embodiments of the disclosed architecture has been evaluated when applied to compensate for the distortion due to a traveling wave tube amplifier (TWTA) and one of its linearized versions. In the discussion that follows, the performance is presented when the TWTA is described in terms of the Saleh model. Similar results were obtained for the case of an otherwise linearized TWTA. For the Saleh model, the amplitude to phase transfer characteristics function φ(r) in equation (4) (with r denoting the normalized, and possibly modified due to the AM-AM compensation, input signal real envelope) may be represented by:

φ(r)=α_(φ) r ²/(1+β_(φ) r ²); α_(φ)=2.5293, β_(φ)=2.8168   (48)

FIG. 10 plots the phase transfer characteristics function φ(r) versus r for the Saleh model. FIG. 10 also shows the result of approximating the function φ(r) by a polynomial g(r) (see equation (12)) of degree 5. The polynomial approximation is given by:

g(r)=−0.6523 r ⁵+3.015 r ⁴−5.1079 r ³+3.4528 r ²−0.0454 r−0.0004   (49)

As may be inferred from FIG. 10, the two graphs are very close with the rms approximation error equal to 0.0012 rad.

FIGS. 11 and 12 plot, as examples, the first and second derivative of the function g(r) denoted by G₁(r) and G₂(r), respectively, and given by equations (13a) and (13b). As described by equations (22)-(24), the function φ(r) may be approximated by a lower order polynomial about some nominal operating point A₀. FIG. 13 plots the result of such an approximation when A₀=1 and a degree 3 polynomial approximation is used. FIG. 13 also plots the error incurred over the interval [−0.5 1.5] with such an approximation. As may be inferred from FIG. 13, the error is relatively small over the interval. The degree 3 polynomial approximation from equations (22)-(24) may be represented by:

p(r−A ₀)=0.4296(r−A ₀)³−0.3035(r−A ₀)²+0.3351(r−A ₀)+0.6618   (50)

FIG. 14 plots the corresponding result with a degree 2 polynomial given by

p(r−A ₀)=−0.3035(r−A ₀)²+0.3351(r−A ₀)+0.6618   (51)

As may be inferred from FIG. 14, the error of approximation remains relatively small even with degree 2 polynomial. This is desirable as the real time correction is made at the sampling rate f_(s), while the coefficients of the polynomial p(r−A₀) need to be updated only once in the interval of N_(s) samples. In most applications, N_(s) is relatively very large number.

The amplitude transfer characteristics function h(x) in equation (4) for the Saleh model may be represented by:

h(x)=α_(a) x/(1+β_(a) x ²); α_(a)=1.9638, β_(a)=0.9945   (52)

The AM-AM distortion is compensated by a function B_(m)(x) of equation (36) where the polynomial B(x) may be an odd degree polynomial obtained by the iterative optimization described by equations (28)-(35). The parameter A_(p) in equation (36) depends upon the degree of polynomial B(x) selected for optimization. FIG. 15 plots the amplitude transfer characteristics function h(x) given by equation (52). FIG. 15 also shows the result with the distortion compensation when the degree of the polynomial B(x) is 7 with 4 nonzero terms. The polynomial B(x) for this case may be represented by:

B(x)=0.7612x ⁷+0.2532x ⁵+0.0835x ³+1.0292x   (53)

The parameter A_(p) for this case is equal to 0.763.

FIG. 16 shows the result of compensation with a degree 9 odd polynomial with 5 non-zero coefficients. The parameter A_(p) for this case is equal to 0.518 and the polynomial B(x) may be represented by:

B(x)=106.9x ⁹+12.05x ⁷+1.357x ⁵+0.1527x ³+1.0168x   (54)

As may be inferred from FIG. 16, for an operating point near saturation, the compensation provides a significantly higher output power compared to that obtained without any compensation. It is to be noted that even with a relatively large coefficient of degree 9 term, the contribution due to this term is relatively small over the range of x between 0 and A_(p)=0.518 where the polynomial B(x) is computed.

The performance of the compensation system may be evaluated when the input to the amplifier is the sum of a wideband signal and a narrow band signal, wherein both the signals are digitally modulated signals with OQPSK modulation. In the simulation example, both the signals are filtered by square-root raised cosine filters. The WB signal has a rate R_(s1)=100 Msps with the square-root raised cosine filter roll-off factor α₁ equal to 0.3. The narrowband (NB) signal with rate R_(s2)=4.44 Msps has the filter roll-off factor α₂=0.21. The separation between the center frequencies of the two signals Δf_(c) is equal to 68.12 MHz, thus avoiding any overlap among their spectrums. Both the signals have the same peak power spectral densities. The ratio of the power levels of the two signals is 22.5 or 13.52 dB. The performance may be obtained in terms of the probability of bit error for the two signals in the presence of the distortion introduced by the amplifier and the receiver thermal noise. The results may be obtained for various values of the output power back off in dB relative to the maximum possible power level at the amplifier input with a tone signal input. Due to the nature of the band limited digital signals, the minimum achievable output power back off is strictly greater than zero.

FIG. 17 plots the amplifier output back off in dB versus the input signal power level when the signal is the sum of two OQPSK signals described above obtained by simulation. The graph of FIG. 17 may then be used to determine the required input power level for any specified output back off.

FIG. 18 plots the probability of bit error P_(e) for the WB signal versus the effective (E_(b)/N₀) for the wide band signal, denoted by (E_(b)/N₀)_(WB) and defined by:

(E _(b) /N ₀)_(WB)=(E _(b) /N ₀)_(T) −P _(L1); (E _(b) /N ₀)_(T)=10 log₁₀[0.5P _(T) T _(s1) /N ₀]  (55)

where P_(T) is the total power at the receiver input, T_(s1)=1/R_(s1) is the symbol period for the WB signal, and P_(L1) is the power sharing loss for the WB signal defined by:

P _(L1)(dB)=−10 log {R _(s1) /[R _(s1)+10^(ΔP) ⁰ ^(/10) R _(s2)]}  (56)

In equation (56), ΔΠ₀ is the difference between the peak power spectral densities of the NB and WB signals. For the simulation example ΔΠ₀=0 dB and P_(L1)=0.19 dB. The effective (E_(b)/N₀) for the narrow band signal denoted by (E_(b)/N₀)_(NB) is similarly defined by:

(E _(b) /N ₀)_(NB)=(E _(b) /N ₀)_(T) −P _(L1) +ΔP ₀   (57)

The bit energy to noise spectral density ratio (E_(b)/N₀)_(NB) will be equal to the actual (E_(b)/N₀) for the NB signal if the amplifier were linear. For the case of nonlinear amplifier considered here, the two will be different; however, (E_(b)/N₀)_(NB) is still a convenient parameter to be used for comparison. For the case of ΔΠ₀=0 dB, (E_(b)/N₀)_(NB)=(E_(b)/N₀)_(WB) from equations (55) and (57).

FIG. 18 plots the probability of bit error P_(e) for the WB signal versus (E_(b)/N₀)_(WB) for the case of no compensation, with only AM-PM compensation, and with both AM-AM and AM-PM compensation. For reference, FIG. 18 also plots the result for the case of a linear amplifier with no distortion. As may be inferred from FIG. 18, there is a drastic performance improvement with the AM-PM compensation, with the performance being close to the ideal case of the linear amplifier. For example, for a P_(e)=10⁻², the improvement is about 5 dB. In this case the additional improvement due to the AM-AM compensation is only marginal, as with this level of back off the amplifier behaves linearly in terms of amplitude transfer characteristics.

FIG. 19 plots the corresponding result for the probability of bit error P_(e) for the NB signal versus (E_(b)/N₀)_(NB). In this case there is a significant improvement of about 1.4 dB due to AM-AM compensation over that due to AM-PM compensation alone. The performance when both AM-AM and AM-PM compensation are used is very close to that for the ideal linear amplifier case. As expected, the performance without any compensation is very poor. The degree N_(a) of the odd polynomial B(x) is equal to 9, which results in about 1.1 dB increase in the output power with the same input power level due to the AM-AM compensation, thus in effect reducing the input power level for the same 3 dB output back off. The result for the case of N_(a)=7 provides marginally higher P_(e) compared to that for N_(a)=9.

FIGS. 20 and 21 plot the probability of bit error results for the WB and NB signals, respectively, when the output back off is 2 dB and with N_(a)=7. As may be inferred from FIG. 20, there is only a small fraction of dB degradation compared to the case of an ideal linear amplifier for the WB signal when only the AM-PM compensation is used or when both AM-PM and AM-AM compensations are used. As expected, the performance with no compensation is extremely poor. As shown in FIG. 21, the AM-AM compensation improves the result over that provided by AM-PM compensation by about 0.4 dB for a P_(e) of about 10⁻⁵. For a P_(e) of about 10⁻², the compensated system requires about 1.1 dB more than the ideal case. The AM-AM compensation increases the output power by about 0.85 dB for the same input power level due to the AM-AM compensation.

FIGS. 22 and 23 plot the probability of bit error results for the WB and NB signals, respectively, when the output back off is 1 dB and with N_(a)=7. For the WB signal, the residual degradation after the compensation for the WB signal is relatively small even with 1 dB output back off.

FIGS. 24 and 25 plot the probability of bit error results for the WB and NB signals, respectively, when the output back off is 0.5 dB and with N_(a)=9. For the WB signal, the residual degradation after the compensation for the WB signal is relatively small even with 0.5 dB output back off. For example, for a P_(e) of 10⁻², the residual degradation in terms of the required (E_(b)/N₀) for the WB signal is only 0.45 dB compared to the ideal linear amplifier. As may be inferred from FIG. 17, the minimum possible output back off for the TWTA without any compensation is 0.35 dB for the type of input signal used in the simulations. It will thus be appreciated that embodiments of the disclosed architecture result in insignificant degradation in performance compared to the ideal linear amplifier even when operating the amplifier only 0.15 dB back off from the maximum possible output level. Reducing the back off by an additional 0.15 dB is not expected to result in any significant performance difference over that for the 0.5 dB back off. The AM-AM compensation provides an increase of 0.22 dB power compared to the case of no AM-AM compensation for the same input power level in both cases. However, the degradation for the NB signal can be significant for 0.15 dB effective back off unless code requiring a channel error rate of about 0.1 is used, which then requires a (E_(b)/N₀)_(NB) of about 2.5 dB.

The performance of the NB signal can be improved by increasing ΔΠ₀ above 0 dB. Such an increase results in a penalty in terms of an effective increase in the required overall (E_(b)/N₀). The overall increase in the required (E_(b)/N₀) due to a higher (E_(b)/N₀) for the NB signal compared to the WB signal by ΔΠ₀ dB may be represented by:

Δ(E _(b) /N ₀)=10 log {1+[(10^(0.1ΔP) ⁰ −1)R _(s2) ]/[R _(s1) +R _(s2)]}  (58)

From equation (58) for the case considered in the simulation example, Δ(E_(b)/N₀) is about 0.19 dB with (E_(b)/N₀)_(NB)−(E_(b)/N₀)_(WB)=3 dB. With the difference in the (E_(b)/N₀) of the two signals equal to 6 dB, Δ(E_(b)/N₀)=0.53 dB. Thus, the amplifier may be operated at saturation even with two signals at its input with only a minimal increase in the overall required (E_(b)/N₀).

It will be appreciated by one of ordinary skill in the art that at least some of the embodiments described herein or parts thereof may be implemented using hardware, firmware and/or software. The firmware and software may be implemented using any suitable computing device(s). FIG. 26 shows an example of a computing device 2600 according to one embodiment. For the sake of clarity, the computing device 2600 is illustrated and described here in the context of a single computing device. However, it is to be appreciated and understood that any number of suitably configured computing devices 2600 can be used to implement a described embodiment. It will also be appreciated that one such device or multiple devices may be shared in a time division multiplex mode among compensators for multiple power amplifiers, as may be the case, for example, in a base station of a mobile communication network. For example, in at least some implementations, multiple communicatively linked computing devices 2600 are used. One or more of these devices can be communicatively linked in any suitable way such as via one or more networks. One or more networks can include, without limitation: the Internet, one or more local area networks (LANs), one or more wide area networks (WANs), or any combination thereof.

In this example, the computing device 2600 comprises one or more processor circuits or processing units 2602, one or more memory circuits and/or storage circuit component(s) 2604 and one or more input/output (I/O) circuit devices 2606. Additionally, the computing device 2600 comprises a bus 2608 that allows the various circuit components and devices to communicate with one another. The bus 2608 represents one or more of any of several types of bus structures, including a memory bus or memory controller, a peripheral bus, an accelerated graphics port, and a processor or local bus using any of a variety of bus architectures. The bus 2608 may comprise wired and/or wireless buses.

The processing unit 2602 may be responsible for executing various software programs such as system programs, applications programs, and/or program modules/blocks to provide computing and processing operations for the computing device 2600. The processing unit 2602 may be responsible for performing various voice and data communications operations for the computing device 2600 such as transmitting and receiving voice and data information over one or more wired or wireless communications channels. Although the processing unit 2602 of the computing device 2600 is shown in the context of a single processor architecture, it may be appreciated that the computing device 2600 may use any suitable processor architecture and/or any suitable number of processors in accordance with the described embodiments. In one embodiment, the processing unit 2602 may be implemented using a single integrated processor.

The processing unit 2602 may be implemented as a host central processing unit (CPU) using any suitable processor circuit or logic device (circuit), such as a general purpose processor. The processing unit 2602 also may be implemented as a chip multiprocessor (CMP), dedicated processor, embedded processor, media processor, input/output (I/O) processor, co-processor, microprocessor, controller, microcontroller, application specific integrated circuit (ASIC), field programmable gate array (FPGA), programmable logic device (PLD), or other processing device in accordance with the described embodiments.

As shown, the processing unit 2602 may be coupled to the memory and/or storage component(s) 2604 through the bus 2608. The bus 2608 may comprise any suitable interface and/or bus architecture for allowing the processing unit 2602 to access the memory and/or storage component(s) 2604. Although the memory and/or storage component(s) 2604 may be shown as being separate from the processing unit 2602 for purposes of illustration, it is worthy to note that in various embodiments some portion or the entire memory and/or storage component(s) 2604 may be included on the same integrated circuit as the processing unit 2602. Alternatively, some portion or the entire memory and/or storage component(s) 2604 may be disposed on an integrated circuit or other medium (e.g., hard disk drive) external to the integrated circuit of the processing unit 2602. In various embodiments, the computing device 2600 may comprise an expansion slot to support a multimedia and/or memory card, for example.

The memory and/or storage component(s) 2604 represent one or more computer-readable media. The memory and/or storage component(s) 2604 may be implemented using any computer-readable media capable of storing data such as volatile or non-volatile memory, removable or non-removable memory, erasable or non-erasable memory, writeable or re-writeable memory, and so forth. The memory and/or storage component(s) 2604 may comprise volatile media (e.g., random access memory (RAM)) and/or nonvolatile media (e.g., read only memory (ROM), Flash memory, optical disks, magnetic disks, and the like). The memory and/or storage component(s) 2604 may comprise fixed media (e.g., RAM, ROM, a fixed hard drive, etc.) as well as removable media (e.g., a Flash memory drive, a removable hard drive, an optical disk). Examples of computer-readable storage media may include, without limitation, RAM, dynamic RAM (DRAM), Double-Data-Rate DRAM (DDRAM), synchronous DRAM (SDRAM), static RAM (SRAM), read-only memory (ROM), programmable ROM (PROM), erasable programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM), flash memory (e.g., NOR or NAND flash memory), content addressable memory (CAM), polymer memory (e.g., ferroelectric polymer memory), phase-change memory, ovonic memory, ferroelectric memory, silicon-oxide-nitride-oxide-silicon (SONOS) memory, magnetic or optical cards, or any other type of media suitable for storing information.

The one or more I/O devices 2606 allow a user to enter commands and information to the computing device 2600, and also allow information to be presented to the user and/or other components or devices. Examples of input devices include data ports, ADCs, DACs, a keyboard, a cursor control device (e.g., a mouse), a microphone, a scanner, and the like. Examples of output devices include data ports, ADCs, DACs, a display device (e.g., a monitor or projector, speakers, a printer, a network card). The computing device 2600 may comprise an alphanumeric keypad coupled to the processing unit 2602. The keypad may comprise, for example, a QWERTY key layout and an integrated number dial pad. The computing device 2600 may comprise a display coupled to the processing unit 2602. The display may comprise any suitable visual interface for displaying content to a user of the computing device 2600. In one embodiment, for example, the display may be implemented by a liquid crystal display (LCD) such as a touch-sensitive color (e.g., 76-bit color) thin-film transistor (TFT) LCD screen. The touch-sensitive LCD may be used with a stylus and/or a handwriting recognizer program.

The processing unit 2602 may be arranged to provide processing or computing resources to the computing device 2600. For example, the processing unit 2602 may be responsible for executing various software programs including system programs such as operating system (OS) and application programs. System programs generally may assist in the running of the computing device 2600 and may be directly responsible for controlling, integrating, and managing the individual hardware components of the computer system. The OS may be implemented, for example, as a Microsoft® Windows OS, Symbian OS™, Embedix OS, Linux OS, Binary Run-time Environment for Wireless (BREW) OS, JavaOS, or other suitable OS in accordance with the described embodiments. The computing device 2600 may comprise other system programs such as device drivers, programming tools, utility programs, software libraries, application programming interfaces (APIs), and so forth.

Various embodiments may be described herein in the general context of computer executable instructions, such as software or program modules/blocks, being executed by a computer. Generally, program modules/blocks include any software element arranged to perform particular operations or implement particular abstract data types. Software can include routines, programs, objects, components, data structures and the like that perform particular tasks or implement particular abstract data types. An implementation of these modules/blocks or components and techniques may be stored on and/or transmitted across some form of computer-readable media. In this regard, computer-readable media can be any available medium or media used to store information and accessible by a computing device. Some embodiments also may be practiced in distributed computing environments where operations are performed by one or more remote processing devices that are linked through a communications network. In a distributed computing environment, program modules/blocks may be located in both local and remote computer storage media including memory storage devices.

Although some embodiments may be illustrated and described as comprising functional component or modules/blocks performing various operations, it can be appreciated that such components or modules/blocks may be implemented by one or more hardware components, software components, and/or combination thereof. The functional components and/or modules/blocks may be implemented, for example, by logic (e.g., instructions, data, and/or code) to be executed by a logic device (e.g., processor). Such logic may be stored internally or externally to a logic device on one or more types of computer-readable storage media. Examples of hardware elements may include processors, microprocessors, circuits, circuit elements (e.g., transistors, resistors, capacitors, inductors, and so forth), integrated circuits, application specific integrated circuits (ASIC), programmable logic devices (PLD), digital signal processors (DSPs), field programmable gate array (FPGA), logic gates, registers, semiconductor devices, chips, microchips, chip sets, and so forth. Examples of software may include software components, programs, applications, computer programs, application programs, system programs, machine programs, operating system software, middleware, firmware, software modules/blocks, routines, subroutines, functions, methods, procedures, software interfaces, application program interfaces (API), instruction sets, computing code, computer code, code segments, computer code segments, words, values, symbols, or any combination thereof. Determining whether an embodiment is implemented using hardware elements and/or software elements may vary in accordance with any number of factors, such as desired computational rate, power levels, heat tolerances, processing cycle budget, input data rates, output data rates, memory resources, data bus speeds, and other design or performance constraints.

It also is to be appreciated that the described embodiments illustrate example implementations, and that the functional components and/or modules/blocks may be implemented in various other ways which are consistent with the described embodiments. Furthermore, the operations performed by such components and/or modules/blocks may be combined and/or separated for a given implementation and may be performed by a greater number or fewer number of components and modules/blocks.

It is worthy to note that any reference to “one embodiment” or “an embodiment” means that a particular feature, structure, or characteristic described in connection with the embodiment is included in at least one embodiment. The appearances of the phrase “in one embodiment” in the specification are not necessarily all referring to the same embodiment.

Unless specifically stated otherwise, it may be appreciated that terms such as “processing,” “computing,” “calculating,” “determining,” or the like, refer to the action and/or processes of a computer or computing system, or similar electronic computing device, that manipulates and/or transforms data represented as physical quantities (e.g., electronic) within registers and/or memories into other data similarly represented as physical quantities within the memories, registers or other such information storage, transmission, or display devices.

While certain features of the embodiments have been illustrated as described above, many modifications, substitutions, changes, and equivalents will now occur to those skilled in the art. It is therefore to be understood that the appended claims are intended to cover all such modifications and changes as fall within the scope of the embodiments.

A first method for pre-processing a signal prior to receipt of the signal by a NLD to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD is disclosed hereinabove. In one embodiment, the first method includes: generating, by at least one processor, a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generating, by the at least one processor, a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generating, by the at least one processor, an AM-AM compensated signal based on the scale factor c_(k); generating, by the at least one processor, a first AM-PM phase compensation value based on the modified amplitude r_(k); and generating, by the at least one processor, an AM-AM/AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value. In certain embodiments of the first method, the signal is a complex baseband signal, and a digitized representation of the signal includes an inphase signal component u_(i,k) and a quadrature signal component v_(i,k). In one such embodiment, the AM-AM compensated signal includes an inphase signal component u_(p,k) and a quadrature signal component v_(p,k), and the step of generating an AM-AM compensated signal includes determining u_(p,k) as the product c_(k)×u_(i,k) and v_(p,k) as the product c_(k)×v_(i,k).

In certain embodiments of the first method, the signal is an RF signal or an IF signal. In one such embodiment, the step of generating an AM-AM compensated signal includes multiplying the signal by an analog representation of the scale factor c_(k).

In certain embodiments of the first method, the step of generating a modified amplitude r_(k) includes evaluating a function B(x_(k)) to generate r_(k). B may be an inverse function of an amplitude transfer function h of the NLD.

In certain embodiments of the first method, the step of generating a scale factor c_(k) includes determining rk/xk.

In certain embodiments of the first method, the step of generating a first AM-PM phase compensation value includes: determining coefficient values of a polynomial φ of degree n_(p) that approximates an amplitude to phase transfer function of the NLD; and evaluating φ(r_(k)), wherein φ(r_(k)) is equal to the first AM-PM phase compensation value. The degree n_(p) of φ may be 2 or 3 in certain embodiments. The step of generating a first AM-PM phase compensation value may further include, in certain embodiments: estimating a nominal operating amplitude A_(0,n) based on an average value of the modified amplitude r_(k) over each of a number of time intervals T_(I); and adaptively updating coefficient values of φ such that φ approximates the amplitude to phase transfer function of the NLD about the nominal operating amplitude A_(0,n).

In certain embodiments, the first method further includes: generating a second AM-PM phase compensation value by evaluating a function {circumflex over (ψ)}_(k)(r_(k)); and generating an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the second AM-PM phase compensation value. {circumflex over (ψ)}_(k) may be a polynomial function of degree n_(q) and model a phase error ψ_(k) caused by at least one of: (1) environmentally-induced variations in an amplitude to phase transfer function of the NLD, (2) time-dependent variations in an amplitude to phase transfer function of the NLD, and (3) an unmodeled nonlinearity of the NLD. Coefficient values of {circumflex over (ψ)}_(k) may be adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k).

In certain embodiments, the first method further includes: iteratively determining coefficient values of B using a gradient algorithm for minimizing error between an input amplitude x and h(B(x)); and terminating the iterative determination when h(B(x)) is substantially linear over a range of x. B may be a polynomial of degree N_(a). In certain embodiments of the first method, the degree N_(a) of B is an odd number. For example, the degree N_(a) of B may selected from the group consisting of 3, 5, 7 and 9 in certain embodiments.

In certain embodiments, the first method further includes limiting B to a constant value A_(p) for values of x>A_(th), wherein A_(th) is a minimum value of x for which h(B(x)) is substantially equal to a saturation amplitude A_(s).

In certain embodiments, the first method further includes, based on changes in an average value of the modified amplitude r_(k) over a number of time intervals T_(I), adaptively modifying the degree N_(a) of B.

A second method for pre-processing a signal prior to receipt of the signal by a NLD to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD is disclosed hereinabove. In one embodiment the second method includes compensating, by at least one processor, an AM-AM distortion of the signal based on a scale factor rk/xk. x_(k) may be an amplitude derived from a digitized envelope R_(i,k) of the signal, and r_(k) may be a modified amplitude determined by a function B(x_(k)). B may be a polynomial function of degree N_(a) and an inverse function of an amplitude transfer function h of the NLD. The second method further includes compensating, by at least one processor, an AM-PM distortion of the signal based on a first AM-PM phase compensation value determined by a function φ(r_(k)). φ may be a polynomial of degree n_(p) approximating the amplitude to phase transfer function of the NLD about a nominal operating amplitude A_(0,n) of r_(k).

In certain embodiments the second method further includes: generating a second AM-PM phase compensation value determined by a function {circumflex over (ψ)}_(k)(r_(k)); and compensating, by at least one processor, the AM-PM distortion of the signal based on the second AM-PM phase compensation value. {circumflex over (ψ)}_(k) may be a polynomial function of degree n_(q) and model a phase error ψ_(k) caused by at least one of (1) environmentally-induced variations in an amplitude to phase transfer function of the NLD, (2) time-dependent variations in an amplitude to phase transfer function of the NLD, and (3) an unmodeled nonlinearity of the NLD. Coefficient values of {circumflex over (ψ)}_(k) may be adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k).

Also disclosed hereinabove is a system that, in one embodiment, includes a non-linear electronic signal processing device for processing an input signal, and a pre-processor in communication with the non-linear electronic signal processing device. The non-linear electronic signal processing device may have a non-linear transfer function. In certain embodiments, the non-linear electronic signal processing device may be a power amplifier. In other embodiments, the non-linear electronic signal processing device may include a cascade of multiple signal processing devices. In one such embodiment, at least one device in the cascade of multiple signal processing devices may be a power amplifier.

The pre-processor may receive a signal, process the received signal, and output the processed signal to an input of the non-linear electronic signal processing device. The pre-processor may include at least one processor programmed to: generate a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generate a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generate an AM-AM compensated signal based on the scale factor c_(k); generate a first AM-PM phase compensation value based on the modified amplitude r_(k); and generate an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value.

In certain embodiments the processor(s) may be further programmed to: store a polynomial function B of degree N_(a) and evaluate B(x_(k)), wherein B(x_(k))=r_(k). B may be an inverse function of an amplitude transfer function h of the non-linear electronic signal processing device. In certain embodiments, the degree N_(a) of B may be an odd number. For example, the degree N_(a) of B may be selected from the group consisting of 3, 5, 7 and 9 in certain embodiments. In certain embodiments, the degree N_(a) of B may be based on changes in an average value of the modified amplitude r_(k) over a number of time intervals T_(I).

In certain embodiments, the processor(s) may be further programmed to: store a polynomial φ of degree n_(p) that approximates an amplitude to phase transfer function of the non-linear electronic signal processing device; and evaluate φ(r_(k)), wherein φ(r_(k)) is equal to the first AM-PM phase compensation value. In certain embodiments, the degree n_(p) of φ may be selected from the group consisting of 2 and 3.

In certain embodiments, the processor(s) may be further programmed to: estimate a nominal operating amplitude A_(0,n) based on an average value of the modified amplitude r_(k) over each of a number of time intervals T_(I); and adaptively update coefficient values of φ such that φ approximates the amplitude to phase transfer function of the non-linear electronic signal processing device about the nominal operating amplitude A_(0,n).

In certain embodiments, the processor(s) may be further programmed to: generate a second AM-PM phase compensation value by evaluating a function {circumflex over (ψ)}_(k) (r_(k)); and generate the AM-AM and AM-PM compensated signal by modifying the phase of the AM-AM compensated signal based on the second AM-PM phase compensation value. {circumflex over (ψ)}_(k) may be a polynomial function of degree n_(q) and model a phase error ψ_(k) caused by at least one of: (1) environmentally-induced variations in an amplitude to phase transfer function of the non-linear electronic signal processing device, (2) time-dependent variations in an amplitude to phase transfer function of the non-linear electronic signal processing device, and (3) an unmodeled nonlinearity of the non-linear electronic signal processing device. Coefficient values of {circumflex over (ψ)}_(k) may be adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k). 

1. A method for pre-processing a signal prior to receipt of the signal by a non-linear device (NLD) to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD, the method comprising: generating, by at least one processor, a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generating, by the at least one processor, a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generating, by the at least one processor, an AM-AM compensated signal based on the scale factor c_(k); generating, by the at least one processor, a first AM-PM phase compensation value based on the modified amplitude r_(k); and generating, by the at least one processor, an AM-AM/AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value.
 2. The method of claim 1, wherein the signal is a complex baseband signal, and wherein a digitized representation of the signal comprises an inphase signal component u_(i,k) and a quadrature signal component v_(i,k).
 3. The method of claim 1, wherein the signal is an RF signal or an IF signal.
 4. The method of claim 1, wherein generating a modified amplitude r_(k) comprises: evaluating a function B(x_(k)) to generate r_(k), wherein B is an inverse function of an amplitude transfer function h of the NLD.
 5. The method of claim 4, further comprising: iteratively determining coefficient values of B using a gradient algorithm for minimizing error between an input amplitude x and h(B(x)), wherein B is a polynomial of degree N_(a); and terminating the iterative determination when h(B(x)) is substantially linear over a range of x.
 6. The method of claim 5, further comprising: limiting B to a constant value A_(p) for values of x>A_(th), wherein A_(th) is a minimum value of x for which h(B(x)) is substantially equal to a saturation amplitude A_(s).
 7. The method of claim 5, wherein the degree N_(a) of B is an odd number.
 8. The method of claim 7, wherein the degree N_(a) of B is selected from the group consisting of 3, 5, 7 and
 9. 9. The method of claim 5, further comprising: based on changes in an average value of the modified amplitude r_(k) over a number of time intervals T_(I), adaptively modifying the degree N_(a) of B.
 10. The method of claim 1, wherein generating a scale factor c_(k) comprises determining rk/xk.
 11. The method of claim 2, wherein the AM-AM compensated signal comprises an inphase signal component u_(p,k) and a quadrature signal component v_(p,k), and wherein generating an AM-AM compensated signal comprises determining u_(p,k) as the product c_(k)×u_(i,k) and v_(p,k) as the product c_(k)×v_(i,k).
 12. The method of claim 3, wherein generating an AM-AM compensated signal comprises multiplying the signal by an analog representation of the scale factor c_(k).
 13. The method of claim 1, wherein generating a first AM-PM phase compensation value comprises: determining coefficient values of a polynomial φ of degree n_(p), the polynomial φ approximating an amplitude to phase transfer function of the NLD; and evaluating φ(r_(k)), wherein φ(r_(k)) is equal to the first AM-PM phase compensation value.
 14. The method of claim 13, wherein generating a first AM-PM phase compensation value comprises: estimating a nominal operating amplitude A_(0,n) based on an average value of the modified amplitude r_(k) over each of a number of time intervals T_(I); and adaptively updating coefficient values of φ such that φ approximates the amplitude to phase transfer function of the NLD about the nominal operating amplitude A_(0,n).
 15. The method of claim 13, wherein the degree n_(p) of φ is selected from the group consisting of 2 and
 3. 16. The method of claim 1, comprising: generating a second AM-PM phase compensation value by evaluating a function {circumflex over (ψ)}_(k)((r_(k)), wherein {circumflex over (ψ)}_(k) is a polynomial function of degree n_(q) and models a phase error ψ_(k) caused by at least one of: (1) environmentally-induced variations in an amplitude to phase transfer function of the NLD, (2) time-dependent variations in an amplitude to phase transfer function of the NLD, and (3) an unmodeled nonlinearity of the NLD, and wherein coefficient values of {circumflex over (ψ)}_(k) are adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k); and generating an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the second AM-PM phase compensation value.
 17. A method for pre-processing a signal prior to receipt of the signal by a non-linear device (NLD) to compensate for AM-AM distortion and AM-PM distortion of the signal by the NLD, the method comprising: compensating, by at least one processor, an AM-AM distortion of the signal based on a scale factor rk/xk, wherein x_(k) is an amplitude derived from a digitized envelope R_(i,k) of the signal, and wherein r_(k) is a modified amplitude determined by a function B(x_(k)), wherein B is a polynomial function of degree N_(a) and an inverse function of an amplitude transfer function h of the NLD; and compensating, by at least one processor, an AM-PM distortion of the signal based on a first AM-PM phase compensation value determined by a function φ(r_(k)), wherein φ is a polynomial of degree n_(p) approximating the amplitude to phase transfer function of the NLD about a nominal operating amplitude A_(0,n) of r_(k).
 18. The method of claim 17, comprising adaptively modifying the degree N_(a) of B based on changes in an average value of the modified amplitude r_(k) over a number of time intervals T_(I).
 19. The method of claim 17, comprising generating a second AM-PM phase compensation value determined by a function {circumflex over (ψ)}_(k)(r_(k)), wherein {circumflex over (ψ)}_(k) is a polynomial function of degree n_(q) and models a phase error ψ_(k) caused by at least one of: (1) environmentally-induced variations in an amplitude to phase transfer function of the NLD, (2) time-dependent variations in an amplitude to phase transfer function of the NLD, and (3) an unmodeled nonlinearity of the NLD, and wherein coefficient values of {circumflex over (ψ)}_(k) are adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k); and compensating, by at least one processor, the AM-PM distortion of the signal based on the second AM-PM phase compensation value.
 20. A system, comprising: a non-linear electronic signal processing device for processing an input signal, wherein the non-linear electronic signal processing device has a non-linear transfer function; and a pre-processor in communication with the non-linear electronic signal processing device, wherein the pre-processor receives a signal, processes the received signal, and outputs the processed signal to an input of the non-linear electronic signal processing device, wherein the pre-processor comprises at least one processor programmed to: generate a modified amplitude r_(k) based on an amplitude x_(k) derived from a digitized envelope R_(i,k) of the signal; generate a scale factor c_(k) based on the modified amplitude r_(k) and the amplitude x_(k); generate an AM-AM compensated signal based on the scale factor c_(k); generate a first AM-PM phase compensation value based on the modified amplitude r_(k); and generate an AM-AM and AM-PM compensated signal by modifying a phase of the AM-AM compensated signal based on the first AM-PM phase compensation value.
 21. The system of claim 20, wherein the least one processor is further programmed to: store a polynomial function B of degree N_(a), wherein B is an inverse function of an amplitude transfer function h of the NLD; and evaluate B(x_(k)), wherein B(x_(k))=r_(k).
 22. The system of claim 21, wherein the degree N_(a) of B is an odd number.
 23. The system of claim 22, wherein the degree N_(a) of B is selected from the group consisting of 3, 5, 7 and
 9. 24. The system of claim 21, wherein the degree N_(a) of B is based on changes in an average value of the modified amplitude r_(k) over a number of time intervals T_(I).
 25. The system of claim 20, wherein the least one processor is further programmed to: store a polynomial φ of degree n_(p), the polynomial φ approximating an amplitude to phase transfer function of the non-linear electronic signal processing device; and evaluate φ(r_(k)), wherein φ(r_(k)) is equal to the first AM-PM phase compensation value.
 26. The system of claim 25, wherein the at least one processor is further programmed to: estimate a nominal operating amplitude A_(0,n) based on an average value of the modified amplitude r_(k) over each of a number of time intervals T_(I); and adaptively update coefficient values of φ such that φ approximates the amplitude to phase transfer function of the NLD about the nominal operating amplitude A_(0,n).
 27. The system of claim 25, wherein the degree n_(p) of φ is selected from the group consisting of 2 and
 3. 28. The system of claim 20, wherein the at least one processor is further programmed to: generate a second AM-PM phase compensation value by evaluating a function {circumflex over (ψ)}_(k)(r_(k)), wherein {circumflex over (ψ)}_(k) is a polynomial function of degree n_(q) and models a phase error ψ_(k) caused by at least one of: (1) environmentally-induced variations in an amplitude to phase transfer function of the non-linear electronic signal processing device, (2) time-dependent variations in an amplitude to phase transfer function of the non-linear electronic signal processing device, and (3) an unmodeled nonlinearity of the non-linear electronic signal processing device, and wherein coefficient values of {circumflex over (ψ)}_(k) are adaptively determined using a gradient algorithm for minimizing error between {circumflex over (ψ)}_(k) and ψ_(k); and generate the AM-AM and AM-PM compensated signal by modifying the phase of the AM-AM compensated signal based on the second AM-PM phase compensation value.
 29. The system of claim 20, wherein the non-linear electronic signal processing device is a power amplifier.
 30. The system of claim 20, wherein the non-linear electronic signal processing device comprises a cascade of multiple signal processing devices.
 31. The system of claim 30, wherein at least one of the cascade of multiple signal processing devices is a power amplifier. 